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Abstract 

We have recently solved the inverse scattering problem for one- 
parameter families of vector fields, and used this result to construct 
the formal solution of the Cauchy problem for a class of integrable non- 
linear partial differential equations in multidimensions, including the 
second heavenly equation of Plebanski and the dispersionless Kadomt- 
sev - Petviashvili (dKP) equation. We showed, in particular, that the 
associated inverse problems can be expressed in terms of nonlinear 
Riemann - Hilbert problems on the real axis. In this paper we make 
use of the nonlinear Riemann - Hilbert problem of dKP, i) to con- 
struct the longtime behaviour of the solutions of its Cauchy problem; 
ii) to characterize a class of implicit solutions; iii) to elucidate the 
spectral mechanism causing the gradient catastrophe of localized so- 
lutions of dKP, at finite time as well as in the longtime regime, and 
the corresponding universal behaviours near breaking. 
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1 Introduction 



It was observed long ago [Ij that the commutation of multidimensional vector 
fields can generate integrable nonlinear partial differential equations (PDEs) 
in arbitrary dimensions. An important example of PDE of this type is the 
dispersionless Kadomtsev-Petviashvili (dKP) equation 

{ut + uu^)^ + Uyy = 0, u = u{x,y,t) eM, x,y,t eR, (1) 

arising from the commutation 

[^1,^2] = (2) 

of the following pair of 1-parameter families of vector fields: 

Li = dy + Xd^ - u^dx, 

L2 = dt + (A2 + u)d^ + {-Xu^ + Uy)dx, 

X being the spectral parameter. 

The dKP equation is the x-dispersionless limit of the celebrated Kadomtsev- 
Petviashvih (KP) equation [2J {ut + uu^ + Uxxx)x + ^Uyy = 0, a = ±1, 
and therefore it describes, for instance, the evolution of small amplitude, 
nearly one-dimensional waves in shallow water |3] near the shore, when the 
x-dispersion can be neglected. It is also a model equation in the description 
of unsteady motion in transonic fiow j4j and in the nonlinear acoustics of 
confined beams [5]. 

We remark that, if also the y - dispersion term Uyy is negligeable in ([1]), 
the dKP equation reduces to the celebrated (1 + 1) - dimensional Hopf equa- 
tion: ut + uux = 0, the integrable prototype model for the description of 
the gradient catastrophe (breaking) of a localized one - dimensional wave. 
Therefore a natural question is whether the dKP equation can be viewed as 
the integrable universal model in the description of the gradient catastrophe 
of a localized two - dimensional wave. As far as we know, no precise results 
are known in this direction, although some formal considerations have been 
made, in the framework of unsteady motions in transonic fiow The fol- 
lowing issues are, in particular, relevant, i) Do localized initial data evolving 
according to dKP break? ii) Do small initial data also break? ii) If break- 
ing occurs, does it take place in a point of the [x, y) - plane or on a line? 
iii) Do the analytic and geometric aspects of the breaking exhibit universal 



2 



features, like in the one- dimensional case? iv) How are the breaking features 
connected with the dKP initial data? In this paper all these questions will 
find a proper answer. 

Dispersionless (or quasi-classical) limits of integrable PDEs, having dKP 
as prototype example, arise in various problems of Mathematical Physics and 
are intensively studied in the recent literature (see, f.i., [7j - [21] )• The Lax 
representation ([3]) and the Hamiltonian formulation for dKP can be found, 
for instance, in [9] and in [10]. An elegant integration scheme for dKP, 
applicable in general to nonlinear PDEs associated with Hamiltonian vector 
fields, was presented in [10] . A nonlinear d - dressing was developed in [15] . 
Special classes of nontrivial solutions were also derived (see, f.i., [H], [TB]). 

The inverse scattering transform for 1-parameter families of multidimen- 
sional vector fields has been developed in [22] (see also [23]). This theory, 
introducing interesting novelties with respect to the classical inverse scatter- 
ing transform for soliton equations [2U |3], has allowed one to construct the 
formal solution of the Cauchy problems for the heavenly equation [25] in [22] 
and for the following novel system of PDEs 



in [26]. The Cauchy problem for the f = reduction of (jlj), the dKP equation 
([T]), was also presented in [26|, while the Cauchy problem for the m = 
reduction of (jl]), an integrable system introduced in [27], was given in [25] . 
Due to the ring property of the space of eigenfunctions associated with vector 
fields, the inverse problem can be formulated in three different distinguished 
ways, one of them being expressed in terms of a linear integral equation. 

In this paper we show that one of these formulations, the nonlinear Rie- 
mann - Hilbert (RH) problem on the real axis, is an efficient way i) to con- 
struct the longtime behaviour of the solutions of the Cauchy problem for 
dKP; ii) to characterize a class of implicit solutions of dKP, parametrized 
by an arbitrary spectral function of a single variable; iii) to establish that 
localized solutions of dKP generically break (if small, they break in the long- 
time regime), and to study in great detail the universal features of the wave 



Vxt + Vyy + UVxx + V:^Vxy - VyVxx = 0, 



(4) 



arising from the commutation of the vector fields 



Li=dy + {X + Vx)dx - u^dx, 

L2 = dt + (A^ + XV:r + U- Vy)dn, + {-Xun, + Uy)dx, 



(5) 
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breaking: the similarity solution before breaking, the vertical inflection at 
breaking, and the development of a three- valued region in the (x, y) - plane 
delimited by a caustic after breaking; iv) to connect the analytic and univer- 
sal features of this gradient catastrophe with the initial data, via the spectral 
transform developed in [26] . 

The paper is organized as follows. In §2 we present the nonlinear RH 
dressing of the system (jlj) and of its dKP reduction ([T]), showing that it 
provides an implicit spectral representation of the solution depending on the 
solution itself, leading to the breaking of localized solutions at flnite t. In 
§3 we make use of such a RH problem to construct the longtime behaviour 
of the solutions of the Cauchy problem for the dKP equation, showing that, 
also asymptotically, the above spectral mechanism for breaking is present. 
In §4 we select a class of spectral data for wich the vector RH problem 
decouples and linearizes, leading to a class of implicit solutions of dKP. In §5 
we study the analytic aspects of the longtime breaking of a localized initial 
condition, showing that this problem is connected to the wave breaking of 
a localized two-dimensional wave evolving according to the 1+1 dimensional 
Hopf equation. In §6 we describe the analytic aspects of the wave breaking 
at flnite time. 



2 Nonlinear RH inverse problem 

From the inverse scattering transform developed in [26j for the system (jlj) 
and for its dKP reduction ([1]), one extracts the following inverse problem. 

Proposition Consider the vector nonlinear RH problem on the real line: 

7f+(A) =7?(7f-(A)), AgM, (6) 

where the solutions vr^(A) = (7rj''(A), 7r^(A)) G are 2-dimensional vector 
functions analytic in the upper and lower halves of the complex A plane, 
normalized in the following way 

nHX)=(^-^'''\-^''~^'''^+diX-'), |A|»1, (7) 

where 

u=]im (A(7rf(A)-A)), (8) 
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and 7?(C) = (7^l(Cl, C2), 7^2(Cl, C2)) G C^, C e are given differentiate 
spectral data. 

Then, assuming that the above RH problem and its linearized form = 
Ja~ are uniquely solvable, where J is the Jacobian matrix of the transforma- 
tion IQ: Jij = d7li/d(j, i,j = 1,2, the solutions 7f^(A) of the RH problem 
dHD are common eigenfunctions of the vector fields (jS]): Ljn^ =0, j = 1,2, 
where 

y, t) = -yu - lim A ( vrf + tn^^ - (x - Ay) ) , (9) 

and u,v are solutions of the nonlinear system (jlj). 

If the spectral data 71(0 satisfy the reahty constraint 



^(7?(C))=C, VCGC2, (10) 

then the solutions u, v are real: m, f G M. If, in addition, the transformation 
( TZiC) is canonical: 

{7^l,7^2}^-'= 1, (11) 

where {■, ■} is the Poisson bracket with respect to the variables (Ci, C2), then 
V = 0, vr^(A) are common eigenfunctions of the vector fields ([3]): Ljtt^ = 
0, j = 1,2 and u satisfies the dKP equation ([1]). 

Proof. One easily shows that the vectors Ljif^, j = 1,2 solve the linearized 
RH problem: (Ljvf"^) = J{LjTf~), and that Ljif^ — as A — > 00. By 
uniqueness, we infer that vf^ are eigenfunctions of the operators Lj: Ljn^ = 
0, j = 1, 2, implying that u and v are solutions of the nonlinear system (jl]). 

From the reality constraint (fTOl) it follows that 7f+ = vf^ , implying m, u G 
R. From the dKP constraint it follows that {vr^^, vr^} = {vrf , n^}, A G M, 
where {■, ■} is now the Poisson bracket with respect to the variables {x. A). 
Since {vrf , vr^} ^ 1, as A — 00, the analyticity properties of the eigenfunc- 
tions imply that {^i, tt^} = {Trjf , 112} = 1- At last, applying Lj, j = 1,2 to 
this equation, on obtains the Hamiltonian constraint v = 0, and system (jlj) 
reduces to the dKP equation. □ 

From the integral equation 

vrr {X) = -XH-Xy + x-2ut + ^J- jr^ho)^^ (^1 ( ^0 , ^2 ( A')) , 

^ ^ \ (12) 
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characterizing the solutions of the RH problem and from the definition 
([HD, one obtains the following spectral characterization of the solution u: 

u = F{x -2ut,y,t) eR, (13) 

where the spectral function F, defined by 

F (e, y,t) = -J §.R, e, t), 7r2-(A; t)) , 

R ^ ^ (,14j 

7^i(Cl,C2) = + ^,(Cl,C2), J = 1,2, 

is connected with the initial data u{x,y,0) via the direct spectral transform 
developed in [26] . 

We remark that, while the inverse spectral transform of most of the known 
integrable PDEs provides a spectral representation of the solution involving, 
as parameters, the space-time coordinates, the inverse problem of dKP pro- 
vides a spectral representation f[T^ . f[T^ of the solution involving, as pa- 
rameter, also the solution u itself, in the combination (x — 2ut). This is the 
spectral mechanism for the breaking of a generic localized initial condition 
at finite time t. We postpone to §5,6 a detailed analysis of the analyticity 
aspects of such a breaking. 

We also remark that the mechanism responsible for this feature is that 
the vector field L2 is quadratic in A and, at the same time, it contains the 
partial differential operator dx. Due to that, the normalization of the ana- 
lytic eigenfunction of L2 involves the coefficient u of the vector field. It is 
easy to see that these properties are shared by the whole dKP hierarchy, as- 
sociated with time operators L„ involving higher powers of A. An example of 
integrable system not associated with the commutation of vector fields, but 
exhibiting an inverse problem in which the solution is implicit, is the Harry 
Dym equation ([291 EO]), connected with the KdV equation by a hodograph 
transformation. 



3 Asymptotics of the dKP equation 

In this section we concentrate on the study of the longtime behaviour of the 
solutions of dKP. Of course, this study is meaningful only if no breaking takes 
place before, at finite time. If the initial condition is, for instance, small, the 
nonlinear term in ([1]) becomes important only in the longtime regime, and 
no breaking takes place before. 
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As we shall see in the following, the spectral mechanism causing the break- 
ing of a localized initial condition evolving according to the dKP equation is 
present also in the longtime regime. 

We study the longtime t » 1 regime in the space regions 

x = x + vit, y = V2t, X — 2ut,vi,V2 = 0(1), f2 7^ 0, t >> 1. (15) 

Then the system of nonlinear integral equations f|T2l) is conveniently 
rewritten in terms of the functions (pj, j = 1,2 defined by 

0i(A) = 7rr(A)+A2t + Ay-(x-2Mt), 02(A) = 7r2-(A) - A, (16) 

as follows 



'('^) ~ 2ni Ir A'-(A-iO) ^ 



R, (^x - 2ut - (A' - A+)(A' - A_)t + 01 (A'), A' + 02 (A')), j 
where 



1,2, 



(17) 



A± = -|±V^^i + f. (18) 

The fast decay, in t, of 0j, j = 1,2, due to the linear growth in t of the 
first argument of Rj, j = I, 2, is partially contrasted if A+ = A_; i.e., on the 
parabola: 



x+^^=x ^v, + ^ = Oj (19) 

of the (x, y) - plane. 

On such parabola the integral equations read 

i?,(x-2Mt- (A' + f)'t + 0i(A'),A' + 02(A')), J = 1,2. ^ ^ 

Since, in this case, the main contribution to the integrals occurs when A' ~ 
—V2/2, we make the change of variable A' = —V2/2 + jj! /y/t, obtaining: 

(21) 
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(22) 



If |A + f2/2| » r^/^ equation implies that 0j(A) = 0{t~^/^), j = 1,2: 

-f + 02(-f + ^)), J = l,2, 
and, via (USD- (HID, that 



u 



-f + 0,(-f + ^) +0,^.. 



(23) 



For A + V2/2 = n/Vi, jJL = 0(1), instead, 0j(A), j = 1,2 are 0(1^ 



-f + 02(-f + ^), J = l,2. 



(24) 



Therefore it is not possible to neglect, in the integral equations fl2TD . 0i_2 
in the arguments of RiXi it follows that these integral equations remain 
nonlinear also in the longtime regime. 

The above results can be summarized as follows. 
On the parabola (fTop . in the space-time regions (fTsp . the longtime behaviour 
of the solution of the dKP equation is given by 

u = -^G (x + ^- 2ut, +0 , (25) 



Vi \ 4:t '2tJ VVt 

where 



G'(^,r/) = -^ j dfxR2(^^ - fx'^ + ai{fx;^,ri),-ri + a2if^;^,ri)^, (26) 



i?2 is the second component of the spectral data R, and ai,2(/^; ^, "'?) are the 
unique solutions of the nonlinear integral equations 



(/^;e,r]) = 2^/^73g^^.(e-/^" + «i(/^';e,r/), 

K ^ (27) 

ri + a2{fJ.';^,ri)j, j = 1,2. 



We remark that the system of integral equations fl27p characterizes the 
solutions of the following vector nonlinear RH problem on the real fi - axis: 



A^i^i;^,v)=( ^~f)+Oifi~'), |/i|»l, 




(28) 



with 




(29) 



Outside the parabola, the solution decays faster than 1 / ^/t. 

We remark that equation ( 125|) has been obtained under the hypothesis 
that both arguments of G are 0(1); then u = 0{1/Vi). It follows that 
2ut = 0{^/i) and, consequently, also x = x + y'^/{4:t) = 0{\/t), to balance 
the term 2ut. 

From equation (l25l) it follows that the spectral mechanism causing the 
breaking of a localized initial condition evolving according to the dKP equa- 
tion is present also in the longtime regime. The analytic aspects of the 
longtime breaking of dKP solutions are illustrated in §5. 

4 Particular solutions of the RH problem 

In this section we construct a class of explicit solutions of the vector nonlinear 
RH problem ([6]) and, correspondingly, a class of implicit solutions of the dKP 
equation parametrized by an arbitrary real spectral function of one variable 
(the general solution would depend on an arbitrary real spectral function of 
two variables). 

Suppose that the components of the spectral data 71 in (Q are given by: 



in terms of the single real spectral function / of a single argument, depending 
on ^1 and (2 only through their product. 
Then the RH problem becomes 



^.(Ci,C2) = 0e^' 



-1, J = 1,2, 



(30) 




A G M 



(31) 
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and the following properties hold. 

i) The reality and the dKP constraints (fTOl) . (fTT|) are satisfied. 

ii) TTj^TT^ = 71^712 ■ Consequently, (vr^vr^) is just a polynomial in A: 

7r+7r+ = TT^TT^ = ~t\^ - + (x - ?>ut)\ - 2yu + 3t d^^Uy = w{X), (32) 

and the vector nonlinear RH problem ([6]) decouples into two scalar, linear 
RH problems: 



(33) 



iii) Since, from (155]) . 

^+e.(-y/+(A) ^ ^-e'(-P/- J = 1, 2, (34) 

where 

also 7r+e^(~)'/+W, j = 1,2 are polynomials in A. We expand them in powers 
of A and introduce the notation 

< A"/ >=^ [ A"/(^(A))rfA, neK (36) 
From the positive power expansion it follows that 



x-2ut-y < f > -t(< Xf > +^ 



(37) 



^+g*/+(«,(A)) _ ^-gi/ {«,(A)) = ^_ < J >^ 

implying the following explicit solution of the Riemann problem: 



^-XH~Xiy + t<f>) + 
X - 2ut - y < f > -t Xf > + 

T,f = {X-<f>) e-*^^(-W). 
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<£>! 

2 



i/±(«,(A)) 



(38) 



From the 1/A terms of 

^+^^i-yf+w^ J = 1,2, that must be zero, one finally 
obtains the two conditions 



u =< Xf > 



2 



2t d~^Uy = yu- {x - 2ui) < f > +y [< Xf > i 



2 




(39) 



t (< AV > + < / X A/ > +^) 



Since, through fl32l) and fl36l) . < A'^/ >, n = 0,1,2 depend implicitely 
on the solution u and on (t d~^Uy), the system ( l39i) characterizes a class of 
implicit solutions of the dKP equation, parametrized by the arbitrary real 
spectral function / of a single variable. 

5 Longtime breaking of dKP solutions 

In this section we show that longtime solutions of dKP break. Of course 
this analysis is meaningful only when breaking does not take place before, 
at finite time. If, for instance, the localized initial condition is small, the 
nonlinear term of the dKP equation becomes relevant only in the longtime 
regime, and breaking can occur only when t is large. 

Let U{x, y, t) be the exact solution of the functional equation (l25il at the 
leading order; i.e.: 



where G, connected with the initial condition of dKP via the direct spectral 
problem introduced in [26], is a largely arbitrary differentiable function of 
two arguments. It is easy to verify that U is the general solution of the first 
order PDE in 2 + 1 dimensions: 



f/(x.,.«) = i=G(x + |!-2t/4). 



(40) 




X 



0. 



(41) 



Equation (17^ suggests to introduce the convenient variables: 



V = ViU, 

ry ryt _l_ ^ 



y = f, i=2Vi, 



(42) 
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transforming the PDE fl4ip into the 1 + 1 dimensional Hopf equation: 



0. 



(43) 



Then the longtime behaviour of the dKP solutions is reduced to the 
study of the evolution of a two-dimensional localized wave under 
the 1 + 1 dimensional Hopf equation ( I43p ; its solution, defined implicitly 
by the equations 

^ = x-G{^,y)t, ^^^> 

depends on y only parametrically. 

Each small portion of the two-dimensional localized wave, characterized 
by its own amplitude V , travels with constant speed V in the x - direction. 
Equations fl44b ) define a 2-parameter family of characteristic curves in the 
(5, t)-space, the parameters being ^, y. Due to the localisation properties 
of G, on each plane y = const, the characteristic curves obtained varying ^ 
intersect, giving rise to the breaking phenomenon. The first breaking will take 
place therefore in a specific point (xf,, yi,, ib) and, going back to the physical 
variables inverting the transformation fl42p . in a specific point {xi,,yb) of the 
{x, y) - plane, at the breaking time tb (see Fig.l). 

5.1 2D - breaking according to the Hopf equation 

In this section we discuss in detail the analytic aspects of the wave breaking 
of a localized two-dimensional wave evolving according to the Hopf equation 
(143|) . whose solution is defined imphcitely by equations (jH]). 

One solves (jHb) with respect to the parameter ^, obtaining C,{x,y,i), 
and replaces it into flUk). to obtain the solution V = G{^{x,y,i),y). The 
inversion of equation fl44b ) is possible iff its ^-derivative is different from 
zero. Therefore the singularity manifold (SM) of the two dimensional Cauchy 
problem for the Hopf equation is the two - dimensional manifold characterized 
by the equation 



S{^,y,t) = l + G^{^,y)i=0 i 



1 



(45) 



Since 




l+Ge^(i,y)i 



(46) 
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the slopes of the locahzed wave become infinity (the so-called gradient catas- 
trophe) on the SM, and the two-dimensional wave "breaks" . 

Then the first breaking time tb, and the corresponding characteristic pa- 
rameters = i^b, Vb) are defined by 

and characterized by the equations: 

Qfe) < 0, G^S) = G^yi^b) = 0, .^g. 

> 0, « ^ G^^S)G^yyi^b) - > 0. ^ ^ 

The corresponding point Xb = {xb, Vb) in which the first wave breaking 
takes place is, from (Hib ): 

Xb = ib + G{(b)ib. (49) 
Now we evaluate equations fl44b) and ( l45l) near breaking, in the regime: 

X = Xb + x\ y = yb + y, i=ib + i' , ^ = ^ + (50) 

where x', y', i', C,' are small. Using fHTj) - fH9l) . equation (jUb) becomes, at the 
leading order, the following cubic equation in ^': 

+ a{y')e' + b{y', i')^ - lX{x', y', i') = 0, (51) 

where 

<y') = 'SS^y'^ biy', i') = ^ [G^e + G^yyy'' , 

X{x', y', i') = S:'- G(a, yb + y')i' - [G(6, yb + y') - G] h ~ (52) 

X' + §|y' - Gi' + ^/ - Gyy'i' + ^y'\ 7 = ?g, 

where 

2t' , , 

6 ^ ^, (53) 

corresponding to the maximal balance 

\a\v\=0{W\ |X| = 0(|e|3/2). (54) 
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In fl52l) and in the rest of this section, all partial derivatives of G whose 
arguments are not indicated are meant to be evaluated at ^f, = {E,b, Vh)- 

The three roots of the cubic are given by the well-known Cardano's for- 
mula: 



e'o {i\y'rt')=-% + {A^)\ + {A.)h^ 



(55) 



where 

A± = R±^/~K (56) 

and the discriminant A reads 

A = + q3^ (57) 

with 

R{x',y',f) = 2x{x',y',i') + ^^ + IQ{y',i'). 
At the same order, function S in fHSj) becomes 

5(e, t) = G^t' + I + 2%se'y' + G^^sy'^l^. (59) 

Known ^' as function of {x,y,t) from the cubic (ISTI) . the solution of 
the Hopf equation and its gradient are then approximated, near breaking, by 
the formulae: 

V{x,y,i)r^G{^, + C,yb + y'), 

^ _Y V(g/,^/)Gfe+g',gft+siO (60) 

Together with the breaking point Xb, another distinguished point is the in- 
flection point Xf-. 

= i^fit')yyb), Xf{i') = Xb + Gi' (61) 
at which R = X = y' = a = C,' = ^'^^ = ^'^y = and, consequently: 

V = G, V(i,^)V^ = i(l,gf), V^^ = V^y = 0. (62) 

We observe that xj = Xb at i = ib- 
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5.1.1 Before breaking 

If t < tb it' < 0), the coefficient Q, defined in f l58|) . is strictly positive, due 
to fj48|) : then the discriminant A = + is also strictly positive and only 
the root ^'q is real. Correspondingly, the real solution of P3l) is single valued 
and described by Cardano's formula. In addition, function S in fISIJl) is also 
strictly positive and V(^x,y)V is finite V x,y. 

To have a more explicit solution, we restrict the asymptotic region around 
the infiection point to the narrow strip : 

|y'| = 0(e«), \X{x',y',i')\=0{e^~^'), (63) 

where 

max (£±i,p) < g < p + 1, p>^, (64) 

obtaining 



X ^x' + {Gy/G^)y' - Gi' = A? ■ (f - f /(?)) = 0{eP+^] 



x'+{Gy/G^)y'-Gi' 

where 



(65) 



/c=(l,gf) (66) 



is normal to the strip and defines the breaking direction. 

Replacing fl^S]) into fpU]) we obtain that, in the narrow strip 0631) - fl64p . 
the solution exhibits a universal behaviour, coinciding with the following 
exact similarity solution of the Hopf equation (H3!) 

Y ^ S:-Xb+iGy/G^)(y-yt) 

t tf) 

and describes the plane tangent to the wave at the infiection point (see Fig.2). 
In addition, the gradient of V behaves like 

V(. y)V , , (68) 

+ 2IG5I I, I' + G««^ ) + slGelG^e^y 

implying that, in the narrow strip fl^ - flM|) . 

V,.,„l'-^^(l.|) (69) 

while outside, in the region |X| = 0(|e|), it becomes finite. We finally remark 
that the tangent plane described by ( 1671) passes through the breaking point. 
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5.1.2 At breaking 

As t t tb, the inflection point becomes the breaking point: Xf — > Xb, the 
above tangent plane (now tangent at the breaking point) becomes vertical, 

with equation k ) = 0, the above strip reduces to the breaking point 

Xb, and the slopes Vj., — > oo as {i — h)'^ (see Fig. 3). 

At the breaking time t = tb, one can give an exphcit description of the 
vertical inflection in the following two subcases. 

i) y' = 0. In this case, the cubic fl5Tl) simplifies to ^'^ = 7x', and the solution 
V exhibits the typical vertical inflection preceeding the wave breaking: 



v.a{,.^VW^,y.) ^v.-lf^^^. (70) 

ii) x' = 0. In this case the cubic simplifies to ^'"^ — '-y{Gy/G^)y' ~ 0, and V 
describes a vertical inflection also wrt y: 



V-G\ib+?h'^{y-yb),yb + y'\ ^ Vy^^-^l^^^. (71) 



G. 




5.1.3 After breaking 

After breaking, the solution becomes three- valued in a domain of the (x, y) 
- plane, and does not describe any physics; nevertheless a detailed study of 
the multivalued region is important, in view of a proper regularization of the 
model, and/or in view of the introduction of a proper single- valued shock 
replacing the multivalued solution. 

lii>ib it' > 0), in the regime ([5lD, the SM equation 5 = 0: 

G^^^e + 2G^^yi'y' + G^yyy'' = \G^\e (72) 

describes an elliptic paraboloid in the (^, y, i) space, with minimum at the 
point (^b,tfc) (see Fig. 4). The intersection of the elliptic paraboloid with any 
t-constant plane, i > ib, defines an ellipse in the {C,,y) - plane (see Fig. 4), 
whose coordinates vary in the following intervals of order 0(|e|^/^): 



1^1 
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Eliminating ^' from equations (1721) and (15T]) . one obtains the SM equation 
in space-time coordinates: 



m\Gh^ (£' + gf - Gi' + ^y'' - Gyy'i' + ^y'') + 

{^e + y'') y' - aC,,, - y'') y'] " = (74) 
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coinciding with the A = condition. It describes a closed curve of the (x, y) 
plane possessing two cusps at the points 



^?(?)~^.tV^^(|.i). (75) 

If {x±{y),y) are the two intersections of this curve with the horizontal 
line y = const, then 

x+(y)-x_(y) = -i|Qp/2 = 0(|e|3/2). (76) 

It follows that this curve is the boundary of a narrow region of thickness 
0(e^/^) in the longitudinal direction, and of thickness 0(e^/^) in the transver- 
sal direction. In Fig. 5 we show its time evolution after breaking. 

Away from the cusps, the solution of the Hopf equation is three valued, 
but two of the branches coincide: 

Vo = G{^, + eo,y), v+ = v^=Gi^, + e+,y) (77) 

and the slopes of the coincident solutions are oo. We remark that the closed 
curve, being the envelope of the intersections with the (x, y) - plane of these 
vertical planes tangent to the wave, is a caustic (see Fig. 6). 

At the two cusps, characterized by the condition Q = R = 0, the three 
real roots of the cubic coincide and the corresponding three solutions of the 
Hopf equation coincide as well: 



Inside the caustic, the discriminant A is strictly negative, the cubic admits 
three different real roots and the solution of the Hopf equation is three- valued. 
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This is the multivalued region that has to be replaced by a proper shock layer, 
whose features depend on the desired regularization. Outside the singularity 
manifold, A > and the solution of the Hopf equation is single valued. 

We end this section remarking that the similarity solution before break- 
ing, the vertical inflection at breaking, and the caustic after breaking make 
clear the universal character of the gradient catastrophe of two-dimensional 
waves evolving according to the Hopf equation. As we shall see in the fol- 
lowing sections, similar considerations hold for the gradient catastrophe of 
two-dimensional waves evolving according to the dKP equation. Similar re- 
sults on a different model can be found in [31]. A general singularity theory 
for caustics and wave fronts can be found in |32j . 



5.2 Longtime breaking of dKP waves 



Inverting the transformation P2l) . the formulas of §5.1 allow one to describe 
the longtime breaking of dKP solutions. Now 



f/(a;,y,t) = ^G(e,y), 



(79) 



and 



Let S,b = {C,b, Vb) be the breaking parameters and (x^, y^, tb) be the break- 
ing point associated with the solution (jUj) of the Hopf equation fH51) . and 
characterized by equations fllTl) - fl^ . Then a localized solution U (x, y, t) of 
equation (!79l) first breaks at the time 



^^=[i) =(2Gcr, (81) 

in the point Xi, = {xb, yt) of the {x,y) - plane given by 

Xb = Xb- vlh, Vb = 2ybtb; (82) 

we remark that Xb is the intersection of the parabola x + y'^/4:tb = Xb with 
the straight line y = yb- 
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The inflection point = {xf{t),yf(t)) is given by 

Xf{t) =Xb + 2G{Vi -Vh)- Vbit - tb) ~ Xb+ 
[2\G^\Gil-e/A)-yl]it-tb), (83) 
yf{t) = yb + 2yb{t - tb), 



(84) 



where the small parameter e, introduced in fl53|) . reads 

t-tb 

e = 

tb 



in terms of the dKP coordinates. The inflection point Xf{t) is the intersection 
of the parabola x + y'^ /At = Xf = Xb + 2G(\/t — y/tb) with the straight line 
y = yf{t) (see Fig. 7), where is deflned in flMl) . At the inflection point, 
Uxx = Uxy = 0, and U and its gradient take the following values: 

U = f,, V(.,,)f/=y|^(l,y, + J,i|)~^(l,y, + 2G,G,). (85) 

Throughout this section we make a systematic use of the expressions of 
the small variables x' , y' ,t' in terms of the dKP variables: 

x' = {x- Xb) + m [{y - yt) - y,{t - tb) + ^) (i - e + 0{e')), 

y, ^ y^ y-y.-2Ut-t,) _ ^ ^ ^(^2)) ^ y,yzMl{l - e + 0{e^)), (§6) 
t'=*^(l-e/4). 

5.2.1 Before breaking 

If t < tf,, in the regime (IMI) . the real solution U of (179!) is single valued and 
described by 

U{x,y,t) ^ ^G{ib + i'o,yb + y) (87) 

where is deflned in f l55l) by Cardano's formula, and x', y' are given in ( l86ll . 

To obtain a more explicit solution, we restrict the analysis to the narrow 
strip around the inflection point, deflned by the conditions 

\X\^\k.{x-Xfm=0{\tr'), ^^ = 0(|e|^), (88) 

where p, q are deflned in flM|) and the vector 

k={l,yb + 2G^Gy) (89) 
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characterizes the breaking direction. On this strip the cubic hnearizes and 
its solution reads 

; _ (x - Xfc) - {2G\G^\ - yl){t - h) + {y, + 2G^Gy){y -y,- 2y,{t - h)) 

^(^-^^) 

(90) 

and U is approximated by the exact similarity solution of equation (HTI) : 

2v^(v^-v^) 
Equation (19T]) is equivalent, at the same order, to 

jj ^ /l^ ^~Xb+vl{t~tb)+{Vb+2G^Gr,)[y-yb-'ivbit-tb)] ^ 

t-tt ' 

describing the plane tangent to the wave at the point Xf (see Fig. 7). The 
gradient of U reads, from ( IHOl) . 



1 (Gg, f^G^ + Yt^y) 



(93) 



where ^' and y' are given in (I90p and fl86p . On the narrow strip (1881) it is 
very large: 

~ {l,yb + 2GyG^) , (94) 

while outside, for |X| = 0(|e|), it is finite. 



5.2.2 At breaking 

As t t tb, the inflection point becomes the breaking point: Xf Xb, the above 
tangent plane becomes vertical, with equation k ■ {x — Xh) and Ux,Uy 00 
as (t — tb)~^. 

At the breaking time t = tb, one can give an explicit description of the 
vertical inflection in the following two subcases. 
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If y = yb, the solution exhibits the typical vertical inflection preceeding 
the wave breaking: 

U ~ -^G Ub + </l{x~Xb), Vb) ^ f/. ~ fl,,,^^ (95) 

li X — Xb + rfbiy — Vb) + = (a line approximately tangent to the 

parabola at Xh), U describes a vertical inflection also wrt y: 

(96) 

5.2.3 After breaking 



If t > tfe, in the regime flMl) . the intersection of the SM with any t - constant 
plane, deflnes an ellipse in the 77) - plane (see Fig. 4), corresponding to the 
following caustic in the (x,?/) - plane, deflned, as in fITIl) . by 

'?,\GAG%.X + ^«^^f ^^^^ (-^e + y'^) y'- 

where now 

X = a; - Xb - [2G'|G'd(l - ^e) - yl{l - e)]{t - tb) + 

[{yb + 2G^Gy){y - yf{t)) + G^y - ybf] (1 - e) + g^y,^ (^)' " 

2|G,|G,^,(t - tbY-^ + ^yl {^y. 

(98) 

The caustic exhibits two cusps at the points 



x^{t) r^XbT \l2^{t - tb). {yb + 2G^Gr„ -1) (99) 
V a 

In addition, if {x^{y),y) are the two intersection points of the caustic with 
the line y = const., we have 

2a3/2 ^\G^\G^^^ fy-y^\\ 3/2 



3 1 1^555 
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therefore the caustic is the boundary of a narrow region of thickness 0(|e|^/^) 

in the longitudinal direction, and of thickness 0(^2^^^(t — 4)) in the transver- 
sal direction (see Fig. 5). On it, the discriminant A of the cubic is zero and, 
away from the cusps, the solution of equation (1791) is three valued, two of the 
branches coincide: 

Uo = ^G{^t + eo,yb + y'), u+ = u^ = ^G{^b + e+,yb + y') (loi) 

and the slopes of the coincident solutions are oo (see Fig. 6). At the two cusps, 
characterized by the condition Q = R = 0, the three real roots of the cubic 
coincide and the corresponding three solutions of equation fl79|) coincide as 
well: 

Uo = U, = U-=^g{^± 2^^^V^(t-t,), ^^^^^ 



Inside the caustic, the discriminant A is strictly negative, the cubic admits 
three different real roots and the solution of equation ( 1791) is three-valued. 
Outside, A > and the solution of the equation fl79l) is single valued. 

The formulae of this section describe, after replacing U by u, the long- 
time breaking of the dKP solutions u if, for instance, the dKP initial data 
Uo{x,y) = u{x,y,0) are small. For small initial data, the inverse spectral 
transform for dKP, developed in [SB], simplifies enormously. The RH spec- 
tral data are expressed in terms of the initial data as follows: 

i?2(Ci,C2) ~ - / ^^uo^i^' + C2y,y), (103) 

m J Ci - f 

and function G, appearing in all formulas of this section, is also given explic- 
itly in terms of uq via fl26|) : 



Gi^^v) ~ ^/ ^^^^^«oe(r- (104) 

Summarizing, we have shown that small and localized initial data evolving 
according to the dKP equation break in the longtime regime, and we have 
described the analytic aspects of such breaking in a surprisingly explicit way. 
The similarity solution before breaking, the vertical inflection at breaking, 
and the caustic after breaking make clear the universal character of such a 
gradient catastrophe. 
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6 The gradient catastrophe at finite time 



Small and localized initial waves evolving according to the dKP equation 
break in the longtime regime in the way described in §5.2. If the initial wave 
is not small, the gradient catastrophe takes place at a finite time and, in this 
section, we study the universal features of this phenomenon immediately 
before breaking. 

The analysis is very similar to that of §5; the main difference is that, while 
function G, appearing in the implicit equation (jUj) for the longtime solution 
of dKP, is an essentially arbitrary smooth function of two variables, function 
F, appearing in the implicit equation f[T51) for the solution of dKP at finite 
time, is a specific function of three variables, connected with the spectral 
data via (fT^ . As we shall see, this implies that, at breaking, F satisfies an 
extra condition (equation fll2ip below) that is instead automatically satisfied 
when F is replaced by its asymptotic form: 

^ i=G(.+ |!,|). (105) 

The inverse problem of the dKP equation, summarized in formulae (fT3!) 
and (fT^ . defines an implicit system of equations for the solution u: 

u = F{C,y,t), 

x = C + 2FiC,y,t)t. ^'^^^ 

In analogy with the considerations of §5, the singularity manifold of dKP 
is the two - dimensional manifold characterized by the equation 

S{C,y,t) = l + 2F^{C,y,t)t = 0, (107) 

that can be solved with respect to t, if + F^tt 7^ 0, giving 

t = i{C,y). (108) 

Since 

the gradient of the localized wave becomes infinity on the SM, and the wave 
"breaks". Let t^, the first time at which 5 = in a point Cfe = {Cb,yb) of the 
(C, y) - plane: 

l + 2Fg(C;,4)t6 = ^ tfe = t(0>); (110) 
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then we obtain the following conditions characterizing the breaking point 
l + 2tfeFf(44) = 

F^iCb,h) < 0, F^jCb,h) + hF^tiCb,h) < 0, ^^^^^ 

FcdCb,^) = -^c yiCb,h) = 0, 

^CCC(0>, h) > 0, P = F^cc(0., tb)Fcyy{Cb, h) - F^^yiC,, u) > 0. 

Due to fll06b ). at the breaking time the wave breaks in the point Xf, = 
{xb, yb) of the (x, y) - plane defined by 

Xb = Cb + 2F{Cb,tb)tb. (112) 

As before, we evaluate equations fll06b ) and (11071) near breaking, in the 
regime: 

X = Xb + x\ y = y^ + y\ t = U + t' , ( = + (' , (113) 
where x',y',t',C ^-^e small, obtaining, at the leading order, the cubic 

C' + a{y')e + h{y\ t')C - lX{x', y', t') = 0, (114) 

where now 

X{x', y', t') =x'- 2F(a, y, t)t' - 2 [F(a, t) - F]tb 

x' + ^y' - 2{F + tbFt)t' - ^y'^ - 2{Fy + tbFy,)y't' - ^^^^^ 

and e is given again by (153|) . corresponding to the maximal balance |C'|, ll/'l = 
0(|e|^/^) and X = 0(|ep/^). In formula (11151) . and in the rest of this section, 
all partial derivatives of F whose arguments are not indicated are meant to 
be evaluated at the breaking point (Cb,4)- 

The three roots of this cubic are given by Cardano's formulas (l55l) -( |571) . 
where now 

Q{y', t') = ^^{F^ + tbFct)e + ^y'^ (116) 
Function S reads, at the leading order, 

S ~ 2(Fc + tbF^t)t' + (FcccC'' + ^F^Cyy'C + ^Cyy?/'') tb- (117) 
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For t < tb {f < 0), since the coefficient Q is strictly positive (see fill 11) 
and (11161) ). so is the discriminant A in (lFr|) . It follows that is real and ^^|_ 
are complex conjugate roots; therefore the solution of dKP is single valued 
and it is described Cardano's formula. In addition, function S is also strictly 
positive, implying that the gradient of u is regular near breaking, in the 
space-time region fl51|l . 

Restricting the asymptotic region in the (x, y) - plane to the narrow strip 
\y'\ = 0{\e\'^), \X\ = 0(|e|P"^^), where p,q are defined by the inequalities 
(IMD, then IC'I = 0{\e\P), C'^.aC,'"^ « bC ~ 7X, and the expressions for X 



;il8) 



and (' simplify: 

X = x' + ^y'- 2{F + UFt)t' + 0{\t\P+^] 



^'+^y'-2{F+tt,Ft)t' 



implying 
Since 



, _ (x - Xb) + {Fy/F^){y - yt,) + c{t - h) 

t — tb 

is an exact similarity solution of dKP, where c is an arbitrary constant, com- 
paring (I119P and (I120p . it follows that F must satisfy the following condition 
at breaking: 

F^ + 2tbF^t = F^t = F^) , (121) 

implying that 

S ~ F^t' + (FcccC' + 2Fayy'C + F^yyy'') h, 

Therefore, in the narrow strip 

\x' + ^y' - 2{F + tbFt)t'\ = 0{\e\P+') (123) 

of the (x, j/)-plane, the solution u is described by the exact similarity solution 
of dKP: 

(x - Xb) + iFy/F^){y - Vb) - (F + 2tfcF,)(t - h) 

U ~ Usim = 7 7 > (-'-24) 

t — tb 
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while the asymptotic expression for the gradient of u follows from equations 
([TO and ^m -. 

^(xv)U 5 . (125) 

pf> I Fccc f x'+iPy/Fc)y'-2iF+ttFt)t' -FcC. V , a „,2 

Then 

V(.,,)n = -^(1,F,/F^) (126) 

in the narrow strip (11231) . while outside, for |X| = 0(|e|), the gradient is fi- 
nite. The vector k = (1, Fy/F(^), orthogonal to the strip, defines the direction 
of breaking. 

In the limit t | tb, the narrow strip ( 1123^ reduces to the breaking point, the 
above tangent plane becomes vertical, with equation F(^{x — Xb)+Fy{y — yb) = 
0, and the slopes Ux, Uy ^ oo at that point as (t — tb)~^ (see Fig. 3). As before, 
we consider two sections in which the description is simpler. 

If y' = 0, the solution of dKP is described by the typical vertical inflection 
at the breaking point Xb- 

F{Cb+ \/l{x-Xb),yb,tb) Ux^^-j=£^==. (127) 

3 ^[x-XbY 

If x' = 0, we have again the vertical inflection at Xb- 



ur^ F \Cb- \ h^{y - yb),yb, h] ^ Uy — ^^^7 



\F 



(128) 
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Fig.l The 2-parameter family of characteristic curves describing the evo- 
lution of 2D-waves according to the 1 + 1 dimensional Hopf equation. 




Fig. 2 In a narrow strip of the plane, the similarity solution describes the 
plane tangent to the wave, immediately before breaking. 
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Fig. 5. Four consecutive snapshots describing the time evolution of the three- 
valued region delimited by the caustic, immediately after breaking. The 
center of each picture is the breaking point. 
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Fig. 6. The three-valued solution after breaking, and the caustic delimiting 
the multivalued region. 




Fig. 7. The two parabolas containing the inflection and breaking points, im- 
mediately before breaking, together with the narrow strip of the similarity 
regime. 
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